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A  method  for  determining  the  ocean-bottom  optical  albedo  Rb  from  in-water  upward  and  downward 
irradiance  measurements  at  a  shallow  site  is  presented,  tested,  and  compared  with  a  more  familiar 
approach  that  requires  additional  measurements  at  a  nearby  deep-water  site.  Also  presented  are  two 
new  algorithms  for  estimating  Rh  from  measurements  of  the  downward  irradiance  and  vertically  upward 
radiance.  All  methods  performed  well  in  numerical  situations  at  depths  at  which  the  influence  of  the 
bottom  on  the  light  field  was  significant.  ©  1999  Optical  Society  of  America 
OCIS  codes:  010.4450,  030.5620,  100.3190. 


1.  Introduction 

For  shallow  ocean  waters,  knowledge  of  the  optical 
bottom  albedo  Rb  is  necessary  to  model  the  underwa¬ 
ter1"2  and  above-water3  light  field,  to  enhance  under¬ 
water  object  detection  or  imaging,4  and  to  correct  for 
bottom  effects  in  the  optical  remote  sensing  of  water 
depth5’6  or  inherent  optical  properties  (IOPs).7’8 
Measurements  of  Rb  can  also  help  one  identify  the 
bottom-sediment  composition,6  determine  the  distri¬ 
bution  of  benthic  algal  or  coral  communities,9  and 
detect  objects  embedded  in  the  sea  floor.  Further¬ 
more,  values  of  Rb,  defined  as  the  upward  irradiance 
emerging  from  the  bottom  divided  by  the  downward 
irradiance  into  the  bottom,  can  be  used  as  an  integral 
test  for  attempted  measurements  of  the  bottom  bidi¬ 
rectional  reflectance  function. 

Although  the  value  of  the  irradiance  ratio  R(z) 
equals  Rb  at  the  bottom,  it  is  not  possible  to  measure 
R(z)  right  at  the  bottom,  and  irradiance  measure¬ 
ments  just  above  the  bottom  are  difficult  to  obtain 
because  of  instrument  self-shadow.  An  estimate  of 
Rb  can  be  made  by  extrapolating  to  the  bottom  mea¬ 
surements  of  R(z)  at  several  depths  z  near  the  bot¬ 
tom.  However,  extrapolation  is  generally  unreliable 


When  this  research  was  performed,  the  authors  were  with  the 
Department  of  Mechanical  Engineering,  University  of  Washing¬ 
ton,  Box  352600,  Seattle,  Washington  98195-2600.  R.  A.  Leath¬ 
ers  is  now  with  the  NASA  Earth  System  Science  Office,  Code  SA00, 
Stennis  Space  Center,  Mississippi  39529.  His  e-mail  address  is 
rleather@ssc.nasa.gov. 

Received  1  February  1999. 

0003-6935/99/153199-07$15.00/0 
©  1999  Optical  Society  of  America 


because  profiles  of  R(z)  typically  vary  sharply  with 
depth  close  to  the  bottom.10  In  vitro  Rb  measure¬ 
ments  of  small  bottom  samples  can  be  obtained  with 
the  method  in  Ref.  5  or  of  larger  bottom  samples  with 
a  spectral  radiometer.  However,  these  processes  are 
time-consuming.  Also  the  in  vitro  value  of  Rb  is  not 
necessarily  equal  to  the  in  situ  value,  and  it  is  not 
clear  how  representative  these  samples  are  of  larger 
spatial  regions  of  interest. 

Because  in  situ  estimates  of  Rb  from  light  measure¬ 
ments  close  to  the  bottom  are  subject  to  small-scale 
horizontal  variability  of  the  bottom,  determining  Rb 
from  measurements  farther  from  the  bottom  may  be 
preferable,  thereby  obtaining  horizontally  averaged 
values  that  are  more  appropriate  for  remote-sensing 
applications  and  one-dimensional  radiative-transfer 
modeling.  Rb  can  be  estimated9  from  an  R(z)  mea¬ 
surement  just  below  the  sea  surface  R{ 0+)  together 
with  simultaneous  measurements  in  nearby  deep  wa¬ 
ter  of  R( 0+)  and  the  downward  diffuse  attenuation 
coefficient.  Similarly,  a  qualitative  algorithm  has 
been  proposed6  and  tested11-12  for  bottom  character¬ 
ization  from  remote  radiance  measurements  at  two 
wavelength  bands  simultaneously  over  both  shallow 
and  deep  water.  The  disadvantage  of  both  of  these 
methods,  however,  is  that  they  require  a  deep-water 
site  nearby  that  has  the  same  water  composition, 
illumination,  and  surface  conditions  as  the  shallow- 
water  site. 

A  new  method  of  solving  the  inverse-radiative- 
transfer  problem  for  determining  Rb  is  proposed  in 
Section  2.  Measurements  of  the  upward  and  the 
downward  irradiances  at  one  wavelength  and  at  least 
two  mid-water-column  depths  at  only  one  site  are 
required.  Also  proposed  are  two  related  algorithms 
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for  estimating  Rb  from  measurements  of  the  verti¬ 
cally  upward  radiance  and  downward  irradiance. 
These  algorithms,  like  the  previously  developed  ones, 
require  knowledge  of  the  measurement  distances 
above  the  bottom.  Results  of  specific  numerical 
tests  for  all  the  Rb  estimation  algorithms  are  pre¬ 
sented  in  Section  3,  and  a  discussion  is  given  in  Sec¬ 
tion  4. 

2.  Theory 

A.  Preliminaries 

We  are  interested  in  the  azimuthally  averaged  radi¬ 
ance  L(z,  p.)  that  satisfies  the  integro-differential 
radiative-transfer  equation: 

(h  —  +  c)l(z,  p)  =  b  |  3(z,  p,  \i')L(z,  p')dp',  (1) 

where  b  and  c  are  the  scattering  and  the  beam  atten¬ 
uation  coefficients,  3  is  the  azimuthally  integrated 
scattering  phase  function,  and  p  is  the  direction  co¬ 
sine  with  respect  to  the  downward  depth  z.  All  the 
quantities  in  Eq.  (1)  implicitly  depend  on  wavelength. 
The  downward  and  upward  irradiances  are  given  by 

Ed(z)  =  2tt  |  pL(z,  p)dp, 

Jo 

E„(z)  =  2tt  |  ||x|L(z,  p)dp.  (2) 

The  irradiance  ratio  is 

R(z)  =  Eu(z)/Ed(z),  (3) 

and  R{zb)  =  Rb  for  water  depth  zh.  The  analogous 
radiance-irradiance  ratio  is 


where  the  water  depth  at  the  shallow  site  of  interest 
is  assumed  to  be  known  and  the  subscript  2nd  de¬ 
notes  measurements  taken  at  a  nearby  deep-water 
site  characterized  by  the  same  illumination,  sea- 
surface  conditions,  and  water  IOPs  as  the  site  of 
interest.  In  the  derivation  of  Eq.  (6),  K  is  taken  to  be 
K{z )  =  \Kd(z)  +  k(z)]/2  averaged  over  depth  in  a 
nonspecified  manner,  where  k  is  the  coefficient  of 
attenuation  in  the  upward  direction  of  upwelling  pho¬ 
tons1415  (from  both  in-water  scattering  and  the  bot¬ 
tom);  a  qualitative  and  quantitative  study  of  k  is 
given  in  Refs.  5  and  15.  In  practice,  however,  the 
value  of  K2nd  is  typically  approximated  by  the  verti¬ 
cally  averaged  Kd[z)  at  the  deep-water  site.4  A  re¬ 
arrangement  of  Eq.  (6)  gives  the  algorithm  evaluated 
by  Maritorena  et  al.5  for  determining  Rb  from  irradi¬ 
ance  measurements  at  two  sites9-16: 

Rb  =  Rindiz)  +  [R(z)  -  R2nd(z)]exp[2(zb  -  z)K2nd\. 

(7) 

Equation  (7)  is  often  written  for  z  =  0+  in  hope  of 
applying  it  to  remote-sensing  applications.  How¬ 
ever,  it  is  valid  for  any  depth  z  and  is  more  accurate 
at  mid-water  depths  than  at  the  surface.  A  difficulty 
with  implementing  Eq.  (7)  is  that  a  deep-water  (2nd) 
site  may  not  be  available  that  matches  the  water, 
surface,  and  illumination  conditions  of  the  shallow- 
water  site. 

An  alternative  shallow-water  model  was  previ¬ 
ously  derived17  from  the  radiative-transfer  equation 
with  the  eigenfunction  expansion  method.  In  this 
derivation  the  light  field  is  approximated  by  the  sum 
of  two  eigenmodes  that  decrease  in  magnitude  with 
distance  away  from  the  surface  and  bottom.  Ex¬ 
pressed  in  the  form  of  Eq.  (7),  this  model  is 

Rb  =  irk  +  \R(z)  -  RJexp[2(z6  -  z)Rk],  (8) 


RL(z)  =  nLu(z)/Ed(z),  (4) 

where  the  vertically  upward  radiance  Lu(z)  =  L(z, 
—  1).  The  factor  ir  is  included  in  the  definition  of 
Rl(z)  so  that  RL(zb )  »  Rb.  For  a  Lambertian  bot¬ 
tom,  RL(zb )  =  Rb. 

With  increasing  depth  in  optically  deep,  spatially 
uniform,  source-free  waters,  R(z),  RL(z),  and  the 
downward  diffuse  attenuation  coefficient, 


Kd(z) 


1  d Ed(z) 
Ed(z)  dz 


(5) 


asymptotically  approach  values  R  ,_,  R  rL,  and  K,_,  re¬ 
spectively.  These  asymptotic  values  are  IOPs  of  the 
water  that  can  be  uniquely  computed  from  b,  c,  and 

3. 13 


B.  Irradiance-Ratio  Approach 
A  well-known  model  for  the  irradiance  ratio  is5-6 


R(z)  =  R2nd(z)  +  [Rb  -  R2nd(z)]exp[-2(z6  -  z)K2nd\, 

(6) 


which  can  be  used  to  estimate  Rb  from  measurements 
at  a  single  site  provided  that  zb  is  known  and  R  ,  and 
K,  can  be  determined.  Note  that  the  asymptotic 
in-water  irradiance  ratio  R,  in  Eq.  (8)  does  not  rep¬ 
resent  the  same  quantity  denoted  by  that  symbol  in 
Refs.  4  and  9,  where  R «,  is  our  R2lld(0+).  Equations 
(7)  and  (8)  were  derived  with  entirely  independent 
approaches.  However,  their  final  forms  are  very 
similar,  and  the  two  theoretically  converge  when 
R2nd  and  K2nd  are  measured  at  large  depths  in  ho¬ 
mogenous,  source-free  waters,  where  R2nd(z)  =  Rrj 
and  K2nd(z)  =  K,  . 

Equation  (8)  provides  a  new  interpretation  of  the 
attenuation  coefficient  in  Eq.  (7).  The  derivation  of 
Eq.  (7)  suggests  that  K  =  (Kd  +  k)/2,  whereas  the 
derivation  of  Eq.  (8)  suggests  that  K  K,_.  Because 
k  >  Kd,4’15  and  therefore  K  >  Kd,  there  is  a  question 
as  to  the  appropriateness  of  taking  K  to  be  the  ver¬ 
tically  averaged  Kd(z).  Because  typically  k  >  K,,_  > 
Kd,  it  follows  that  K  =  \Kd(z)  +  k(z)]2„(//2  «=  K,_. 
Therefore  Eq.  (7)  may  be  best  implemented  by  taking 
K  equal  to  the  value  of  Kd(z)  deep  in  the  euphotic 
zone,  where  Kd(z)  «<  K,_  rather  than  as  earlier  pro- 
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posed.9"16  This  hypothesis  is  further  addressed  in 
Section  3. 

If  R  ,  and  Kv  can  be  determined  at  the  shallow  site 
of  interest,  the  method  of  Eq.  (8)  has  the  distinct 
advantage  over  that  of  Eq.  (7)  that  measurements  are 
required  from  only  one  site.  One  way  to  determine 
R  ,_  and  K,_  is  to  calculate  them  with  the  procedure  in 
Ref.  17  from  measurements  of  b  and  c  obtained,  for 
example,  from  water  samples18  or  with  a  Wetlabs 
ac-9  instrument.19  Alternatively,  it  is  possible  to  es¬ 
timate  R  ,  and  K,  from  the  same  irradiance-profile 
measurements  (at  the  shallow-water  site)  used  to 
form  R(z)  without  any  direct  measurement  of  the 
water  properties.  The  value  of  Rx  can  be  obtained 
with  an  equation  derived17  from  a  shallow-water  as¬ 
ymptotic  approximation  of  the  light  field: 

ll-  R„\2  =  [Ed(z)-EJz)]%\ 

\1  +Rj  [EJz)  +  EJz)]2 

To  employ  Eq.  (9),  one  must  subtract  and  add  EJz) 
and  EJz)  at  two  depths,  z1  and  z2,  square  the  results, 
and  evaluate  the  differences  between  the  two  depths. 
Because  these  operations  are  susceptible  to  noise,  it 
is  important  that  the  irradiance  measurements  be  of 
high  quality  and  their  temporal  variations  be  aver¬ 
aged  out.  The  value  of  Kx  in  Eq.  (8)  can  be  estimated 
as  the  maximum  value  attained  by  Kd{z).  Because 
the  value  of  Kd(z)  is  relatively  insensitive  to  Rh,]  the 
value  of  ma x.[KJz)]  is  typically  approximately  equal 
to  Kx. 

C.  Radiance-lrradiance  Ratio  Approach 

If  Lu(z)  measurements  are  available  rather  than 
EJz),  Rb  can  be  estimated  from  RL(z)  of  Eq.  (4)  with 
a  new  model  (derived  in  Appendix  A)  analogous  to  Eq. 
(8): 

Rb  =  R,l  +  [RL(z)  -  R./J]exp[2(z6  -  z)KJ\.  (10) 

However,  because  an  equation  analogous  to  Eq.  (9) 
has  not  been  derived  for  R,h,  implementation  of  Eq. 

(10)  requires  that  the  value  of  R  jJ  either  be  calculat¬ 
ed13  from  local  measurements  of  the  IOP’s  or  be  mea¬ 
sured  in  nearby  deep  water. 

We  found  with  numerical  simulations  that  reason¬ 
able  estimates  of  Rb  can  alternatively  be  obtained 
from  RL(z)  with 

Rb  =  R-indiz)  +  \RL{z)  ~  i?2«/(z)]ex p[2(zfc  -  z)K2nd], 

(11) 

although  we  have  no  analytical  justification  for  Eq. 

(11)  other  than  its  analogy  to  Eq.  (7).  As  with  Eq.  (7) 
use  of  this  equation  requires  measurements  at  a  sec¬ 
ond  (deep-water)  site. 

3.  Numerical  Tests 

A.  Methods 

Numerical  tests  were  performed  to  evaluate  the  ac¬ 
curacies  of  Eqs.  (7) — (11)  for  determining  Rh.  Simu¬ 
lated  Eu{z),  Ed{z),  and Lu{z)  values  were  generated  at 


0.25  optical  depth  spacing  by  using  the  discrete  ordi¬ 
nates  radiative-transfer  code  disort.20  The  surface 
illumination  was  modeled  as  a  combination  of  direct 
collimated  sunlight  and  diffuse  skylight.  The  water 
was  defined  to  have  locally  homogenous  optical  prop¬ 
erties,  a  relative  index  of  refraction  of  1.34  with  re¬ 
spect  to  air,  and  scattering  determined  by  the  Petzold 
particle-scattering  phase  function.21  Spatially  de¬ 
pendent  internal  sources,  such  as  from  fluorescence, 
Raman  scattering,  or  bioluminescence,  were  ne¬ 
glected.  A  Lambertian  bottom  was  assumed,  which 
provides  a  good  approximation  to  the  more  general, 
but  usually  poorly  known,  bidirectional  reflectance 
function.1"2  Simulations  were  performed  for  various 
values  of  the  single-scattering  albedo  (w0  =  b/c),  Rb, 
percent  direct  sunlight,  and  water  optical  depth  (t6  = 
czb).  The  values  of  Rb  were  taken  to  be  in  the  range 
0  <  Rb  <  0.4,  which  is  consistent  with  observations  in 
natural  waters.1 

Equation  (8)  was  applied  to  shallow-water  simula¬ 
tions  of  EJz)  and  EJz)  to  determine  bottom  albedo 
estimates  RJz) as  a  function  of  the  depth  of  the  cor¬ 
responding  irradiance  measurements.  The  values 
of  Rr_  and  K,  used  in  Eq.  (8)  were  obtained  in  two 
different  ways.  First,  R  ,_  was  determined  with  Eq. 
(9)  and  K,_  was  estimated  from  either  KJz)  «=  KJz)  or 
Kx  ma -x\KJz)].  Second,  the  IOP’s  of  the  water 
were  assumed  to  be  known  from  in  situ  measure¬ 
ments,  and  these  IOP’s  were  used  to  calculate  R  ,_  and 
K,,  For  comparison,  RJz)  was  also  determined 
from  EJz)  and  EJz)  with  Eq.  (7)  for  combinations  of 
shallow  and  deep-water  simulations.  Here  K  was 
taken  to  be,  alternatively,  KJz),  KJz)  [the  vertical 
average  of  KJz)  between  the  surface  and  the  depth  of 
the  shallow-water  site],  and  Kf  [determined  from  a 
deep  (asymptotic)  value  of  KJz)]. 

The  bottom  albedo  was  also  determined  from  sim¬ 
ulations  of  LJz)  and  EJz).  Equation  (10)  was 
used  to  calculate  RJz)  from  data  at  only  the 
shallow -water  site.  The  values  of  Rj'  and  K „  were 
calculated  from  the  known  water  optical  properties. 
In  addition  Rb(z)  was  determined  from  LJz)  and 
EJz)  with  Eq.  (11)  for  combinations  of  shallow-  and 
deep-water  sites.  The  value  of  K2nd  was  deter¬ 
mined  in  the  same  manner  as  for  the  Eu  —  Ed 
approach. 

B.  Results 

For  all  the  shallow-water  simulations  performed,  es¬ 
timates  of  Rb  obtained  with  Eqs.  (8)  and  (9)  ap¬ 
proached  the  correct  value  in  a  nearly  linear  fashion 
within  the  bottom  few  optical  depths  of  the  water 
column.  Extrapolation  of  RJz)  at  two  and  one  op¬ 
tical  depths  above  the  bottom  consistently  produced 
estimates  of  Rb  thaHwere  accurate  to  within  ~1%. 
In  Table  1  example  RJz)  is  given  at  three,  two,  and 
one  optical  depths  above  the  bottom  as  well  as  the 
linear  extrapolation  of  the  latter  two  to  the  bottom. 
For  Table  1  the  illumination  conditions  of  the  simu¬ 
lations  were  taken  to  be  either  overcast  (100%  dif¬ 
fuse)  or  sunny  (75%  direct  sunlight  from  a  zenith 
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Table  1.  Calculations  of  the  Bottom  Albedo° 


Simulation 

7^(t) 

O)0 

Tft 

R„ 

Illumination 

-  3 

U  ~  2 

T*  -  1 

Extrapolation 

0.7 

3 

0.2 

Sunny 

N/A 

0.177 

0.189 

0.200 

0.7 

5 

0.1 

Sunny 

0.091 

0.091 

0.096 

0.101 

0.7 

5 

0.2 

Sunny 

0.170 

0.176 

0.189 

0.202 

0.7 

5 

0.2 

Overcast 

0.172 

0.177 

0.190 

0.202 

0.7 

7 

0.2 

Sunny 

0.167 

0.176 

0.189 

0.203 

0.8 

5 

0.2 

Sunny 

0.186 

0.188 

0.194 

0.201 

0.9 

5 

0.1 

Sunny 

0.114 

0.106 

0.102 

0.098 

0.9 

5 

0.2 

Sunny 

0.207 

0.201 

0.200 

0.199 

0.9 

5 

0.2 

Overcast 

0.207 

0.201 

0.200 

0.199 

0.9 

5 

0.3 

Sunny 

0.296 

0.294 

0.297 

0.300 

0.9 

5 

0.4 

Sunny 

0.385 

0.387 

0.395 

0.403 

“Calculations  are  with  Eqs.  (8)  and  (9)  from  simulated  irradiance 
measurements  at  three,  two,  and  one  optical  depths  above  the 
bottom  and  extrapolations  to  the  bottom  of  the  latter  two. 


angle  of  30°),  and  in  the  solutions  for  Rb  the  values  of 
K,_  in  Eq.  (8)  were  approximated  by  max[ffr/(z)J.  At 
a  given  measurement  depth  estimates  generally  im¬ 
proved  with  the  increasing  value  of  w0.  For  exam¬ 
ple,  for  Rb  =  0.2,  t6  =  5,  and  sunny  skies,  the  error  at 
three  optical  depths  above  the  bottom  was  15%  for  w0 
=  0.7  but  only  3.5%  for  w0  =  0.9.  The  value  otRb  had 
a  small  effect  on  the  accuracy  of  its  estimate.  For 
large  co0  and  small  Rb,  Rb(z)  was  greater  than  Rb, 
whereas  for  small  values_of  w0  or  for  a  combination  of 
large  w0  and  large  Rb,  Rb(z)  was  less  than  Rb.  Es¬ 
timates  were  more  accurate,  but  insignificantly  so, 
with  overcast  conditions.  Also,  the  depth  of  the  wa¬ 
ter  had  little  effect  on  the  accuracy  ofi?6(z)  at  a  given 
depth  above  the  bottom.  However,  in  practice,  in¬ 
strument  noise  will  be  more  significant  in  relatively 
deep  water  than  in  shallow  water. 

In  most  cases,  estimates  of  Rb  with  Eqs.  (8)  and  (9) 
were  more  accurate  if  K,  in  Eq.  (8)  was  approximated 
by  max[Ad(z)]  than  if  it  was  replaced  by  Kd(z),  be¬ 
cause  large  values  of  K,  in  Eq.  (8)  lead  to  large  values 
of  Rb,  and  the  estimates  oiRb  were  typically  less  than 
the  correct  value.  For  cases  in  which  the  value  of  w0 
was  high  and  the  value  of  Rb  small,  the  use  of  Kd{z) 
gave  slightly,  but  insignificantly,  better  results  than 
the  use  of  max[ffd(z)]. 

Because  Eq.  (9)  provides  only  an  asymptotic  ap¬ 
proximation  to  Rv,  it  was  expected  that  estimates  of 
Rb  would  improve  if  more  accurate  values  of  Rm  cal¬ 
culated  from  the  assumed  known  water  IOP’s,  were 
used  in  Eq.  (8).  However,  the  numerical  tests 
showed  the  reverse  to  be  true;  errors  in  Rb{z)  intro¬ 
duced  by  the  approximation  of  Eq.  (9)  helped  coun¬ 
teract  errors  in  Rb(z)  due  to  the  assumptions 
inherent  in  Eq.  (8).  Although  Eq.  (8)  with  a  calcu¬ 
lated  R  ,_  performed  similarly  in  the  bottom  half  of  the 
water  column  to  Eq.  (8)  with  from  Eq.  (9),  Eq.  (8) 
with  a  calculated  R  ,  performed  poorly  in  the  top  half 
of  the  water  column  and  even  near  the  bottom 
slightly  underperformed  Eqs.  (8)  and  (9). 

For  example,  shown  in  Fig.  1  are  the  estimates  of 


Fig.  1.  Profiles  of  estimates  of  the  bottom  albedo  lFb(  t)  with,  a,  the 
estimate  equal  to  the  irradiance  reflectance  R( t);  b,  Eq.  (8)  with 
known  and  K,,;  and  c,  Eqs.  (8)  and  (9).  The  simulation  was 
generated  for  Rb  =  0.2,  five  optical  depth  water,  w0  =  0.8,  and  75% 
collimated  light  at  30°  from  the  zenith. 


Rb  as  a  function  of  the  measurement  optical  depth 
obtained  from  simulated  Eu{z)  and  Ed{z)  with  three 
different  methods:  from  Rb(z)  =  R{z),  from  Eq.  (8) 
with  calculated  R^  and  K,  ,  and  from  Eqs.  (8)  and  (9) 
with  K,_  replaced  by  max[l?d(z)].  In  this  simulation 
Rb  =  0.2,  t6  =  5,  and  w0  =  0.8,  and  the  sea-surface 
illumination  was  taken  to  be  sunny  (as  defined 
above).  Since  Rb  =  R{z  — >  zb),  the  first  approach  for 
determining  Rb  is  the  most  straightforward.  This 
approach  gave  a  smooth  profile  of  Rb( t)  that  mono- 
tonically  approached  Rb  with  increasing  depth. 
However,  this  estimate  is  extremely  inaccurate  ex¬ 
cept  very  close  to  the  bottom,  with  41%  error  at  only 
one  optical  depth  above  the  bottom  and  61%  error  at 
two  optical  depths  above  the  bottom.  This  sharp 
increase  in  R(z)  near  the  bottom  is  typical10  and 
makes  extrapolation  of  R(z')  from  mid-water  depths 
to  the  bottom  impractical.  At  all  depths  off  the  bot¬ 
tom,  far  better  estimates  of  Rb  were  obtained  with  Eq. 
(8)  with  R  ,_  and  Kr  calculated  from  the  known  water 
IOP’s.  The  errors  at  two  and  one  optical  depths  off 
the  bottom  were  11%  and  5.8%,  respectively.  Even 
better  estimates  of  Rb  at  all  depths  off  the  bottom, 
however,  were  obtained  from  Eq.  (8)  with  deter¬ 
mined  with  Eq.  (9)  and  K,_  estimated  by  max[Arf(z)]. 
This  method  gave  an  error  at  two  and  one  optical 
depths  off  the  bottom  of  6.2%  and  2.8%,  respectively. 
Again,  extrapolation  from  mid-depths  gave  an  excel¬ 
lent  estimate  of  Rb. 

Values  of  the  bottom  albedo  from  irradiance  mea¬ 
surements  at  a  combination  of  shallow-  and  deep¬ 
water  sites  with  Eq.  (7)  were  consistently 
underestimated.  Therefore,  since  K,  was  generally 
larger  than  Kd(z),  results  were  more  accurate  when 
the  value  of  K  in  Eq.  (7)  was  determined  from  a  deep 
value  of  Kd{z)  in  the  deep-water  site,  where  Kd(z) 
K,,  than  when  it  was  determined  from  either  the 
deep-water  Kd[z)  or  Kd(z).  For  example,  shown  in 
Fig.  2  are  Rb( t)  calculated  from  Eq.  (7)  with  K  re¬ 
placed  by  Kd(z)  and  by  K,_  [determined  from  Kd( t  = 
15)]  for  two  sunny-sky  simulations:  (1)  co0  =  0.7,  t6 
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R,  estimate 
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Fig.  2.  Profiles  of  estimates  of  the  bottom  albedo  with  the 
two-site  Eu  -  Ed  method  of  Eq.  (7)  (left)  for  u)0  =  0.7  and  Rh  =  0.2 
and  (right)  for  u>0  =  0.9  and  Rb  =  0.1.  The  simulations  were  for 
sunny  conditions,  and  the  Rb  determination  was  done  with  K in  Eq. 
(7)  replaced  ( —  *  — )  by  Kd(z)  and  ( — )  by  the  deep  (asymptotic) 
value  of  Kd(z). 


=  5,  and  Rb  =  0.1  and  (2)  w0  =  0.9,  t6  =  5,  and  Rb  = 
0.2.  _ 

Shown  in  Fig.  3  are  comparisons  of  7?,,(T)  obtained 
from  Eu(z )  and  Ed(z)  measurements  at  only  the 

shallow-water  site  [Eqs.  (8)  and  (9)]  and  !Fb( t)  ob¬ 
tained  from  measurements  at  both  shallow-  and 
deep-water  sites  [Eq.  (7)].  The  two  cases  shown  are 
sunny-sky  simulations  with  ih  =  5  and  (1)  w0  =  0.9 
andi?6  =  0.1  and  (2)  co0  =  0.8  andf?6  =  0.2.  As  these 
examples  demonstrate,  the  two-site  method  was  typ¬ 
ically  far  more  accurate  than  the  one-site  method 
near  the  surface,  whereas  the  one-site  algorithm  usu¬ 
ally  outperformed  the  two-site  method  near  the  bot¬ 
tom.  _ 

Example  Rb( t)  obtained  from  Lu(z)  and  Ed(z)  are 
shown  in  Fig.  4  for  (1)  w0  =  0.9,  Rb  =  0.1,  and  sunny 
conditions  and  (2)  oo0  =  0.9,  Rb  =  0.2,  and  overcast 
conditions.  Because  estimates  of  Rb  with  Eq.  (11) 
were  typically  larger  than  the  actual  value,  Rb{ t)  in 
Fig.  4  obtained  with  Eq.  (11)  was  calculated  with 
Kd(z)  rather  than  the  asymptotic  Kd  to  make  the 


R,  estimate 
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Fig.  3.  Profiles  of  estimates  of  the  bottom  albedo  Rb(j)  from  Eu(z) 
and  Ed(z)  with  ( —  *  — )  the  two-site  method  of  Eq.  (7)  and  ( — )  the 
one-site  method  of  Eq.  (8).  These  sunny-sky  simulations  were 
(left)  for  o)0  =  0.9  and  Rh  =  0.1  and  (right)  for  <o0  =  0.8  and  Rb  = 
0.2. 


R  estimate 

b 

Fig.  4.  Profiles  of  estimates  of  the  bottom  albedo  ITb( t)  from  Lu(z) 
and  Ed{z)  with  ( —  *  — )  the  two-site  method  of  Eq.  (11)  and  ( — )  the 
one-site  method  of  Eq.  (10).  The  simulations  were  (left)  for  sunny 
conditions,  co0  =  0.9  and  Rb  =  0.1,  and  (right)  for  overcast  condi¬ 
tions,  u>0  =  0.9  and  Rb  =  0.2. 


estimated  values  as  good  as  possible.  In  general  the 
value  of  oj0  had  a  significant  effect  on  the  accuracy  of 
both  Lu  -  Ed  methods,  with  the  best  estimates  ob¬ 
tained  when  co0  was  large.  The  illumination  condi¬ 
tions  were  also  very  important  for  the  Lu  —  Ed 
methods,  with  the  best  results  obtained  for  overcast 
conditions.  The  value  of  Rb,  on  the  other  hand,  had 
no  significant  effect  on  the  accuracy  of  its  estimate. 
In  sunny  conditions  the  one-site  method  of  Eq.  (10) 
and  the  two-site  method  of  Eq.  (11)  performed  simi¬ 
larly  in  the  bottom  one  optical  depth  when  co0  was 
small  and  in  the  bottom  three  optical  depths  when  oj0 
was  large,  but  Eq.  (11)  was  considerably  more  accu¬ 
rate  than  Eq.  (10)  near  the  surface.  For  overcast 
conditions  Eq.  (10)  performed  well  at  all  depths  but 
still  underperformed  Eq.  (11). 

4.  Discussion 

Several  methods  were  evaluated  here  for  determin¬ 
ing  Rb  from  common  natural-light  measurements. 
Each  method  returns  accurate  values  of  Rb  if  imple¬ 
mented  close  to  the  bottom.  However,  because  it  is 
difficult  in  practice  to  obtain  light-field  measure¬ 
ments  close  to  the  bottom,  it  is  necessary  to  apply 
these  algorithms  at  one  or  more  optical  depths  off  the 
bottom  and,  when  possible,  extrapolate  the  depth- 
dependent  estimates  to  the  bottom.  Therefore  it  is 
desirable  that  the  error  of  the  method  used  be  both 
small  and  linearly  decreasing  with  depth.  The 
methods  were  all  found  to  be  preferable  to  a  straight¬ 
forward  extrapolation  of  R{z)  to  the  bottom,  where 
R(z)  =  Rb,  but  differed  in  their  accuracies  when  ap¬ 
plied  more  than  one  or  two  optical  depths  away  from 
the  bottom. 

Estimates  of  Rb  can  be  obtained  with  Eq.  (8)  first  by 
estimating  R^  either  with  Eq.  (9)  or  by  calculating  it 
from  known  water  IOP’s.  This  method  does  not  re¬ 
quire  measurements  at  other  wavelengths  or  at  an¬ 
other  site,  and  estimates  of  Rb  from  this  method  at 
one  or  two  optical  depths  off  the  bottom  were  gener¬ 
ally  found  to  be  more  accurate  than  those  obtained 
with  Eq.  (7).  Numerical  simulations  indicated  that 
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the  use  of  R  ,_  from  Eq.  (9)  produces  better  estimates 
of  Rb  than  the  use  of  the  value  of  R„  computed  from 
the  IOP’s,  because  fortuitously  the  error  introduced 
by  applying  Eq.  (8)  where  the  light  field  is  not  well 
described  by  the  assumed  two-mode  asymptotic 
model  (see  Appendix  A)  is  mitigated  by  the  deviation 
in  the  value  of  R  v  predicted  by  Eq.  (9)  from  its  true 
value.  Therefore,  even  if  the  water  IOP’s  are 
known,  employing  Eq.  (9)  is  preferable,  provided  the 
processed  data  have  a  relatively  smooth  Rb(j)  profile. 
Unfortunately  this  method  is  often  inaccurate  near 
the  sea  surface  when  the  bottom  signal  is  not  strong. 

If  one  wishes  to  estimate  Rb  from  measurements 
close  to  the  surface  and  a  suitable  deep-water  site  is 
available,  Eq.  (7)  provides  the  most  reliable  method. 
However,  it  was  found  that  if  the  deep-water  site  is 
vertically  well  mixed,  Eq.  (7)  should  be  implemented 
by  replacing  K2nd  with  Kr  ,  which  can  be  directly  mea¬ 
sured  deep  in  the  euphotic  zone  of  the  deep-water 
site. 

Estimates  of  Rb  alternatively  can  be  made  from 
measurements  of  Lu(z)  and  Ed(z)  with  Eq.  (10),  pro¬ 
vided  that  the  bottom  is  approximately  Lambertian. 
Equation  (10)  requires  that  R,L  either  be  analytically 
computed  from  local  measurements  of  the  water 
IOP’s  or  be  measured  at  a  nearby  deep-water  site.  If 
a  suitable  second  site  is  readily  available,  Eq.  (11) 
should  be  used  instead  since  it  was  found  to  be  gen¬ 
erally  more  accurate  and  reliable  than  Eq.  (TO). 
However,  both  Lu  -  Ed  methods  performed  well 
when  applied  in  the  bottom  half  of  the  water  column. 
The  inaccuracy  of  Eq.  (10)  in  sunny  conditions  makes 
it  unsuitable  for  remote-sensing  applications,  and 
therefore  Eq.  (10)  offers  no  advantage  over  the  Eu  - 
Ed  method  (which  is  the  more  natural  approach  to  an 
in  situ  Rb  estimation).  Equation  (11),  on  the  other 
hand,  shows  some  promise  for  remote-sensing  appli¬ 
cations  for  large  w0  (or  for  any  co0  if  co0  is  known). 

Regardless  of  the  method  used,  the  determination 
of  Rb  is  easiest  when  the  bottom  signal  is  strong. 
Thus  Rb  can  be  obtained  most  accurately  when  the 
water  is  shallow,  the  value  of  Rb  is  large,  and  the 
attenuation  of  the  water  is  low  (e.g.,  over  tropical 
coral  reefs  or  white  sandy  beaches).  If  the  bottom 
composition  is  believed  to  be  uniform  over  a  large 
horizontal  region,  the  determination  of  Rb  should  be 
made  at  the  shallowest  depth. 

In  practice  the  method  to  use  for  estimating  Rb  will 
be  dictated  by  the  instrumentation  available  and 
whether  an  appropriate  deep-water  site  exists. 
Given  the  choice,  however,  the  estimation  of  Rb 
should  be  made  with  measurements  of  Eu(z)  rather 
than  with  Lu(z).  The  most  informative  approach 
would  be  to  use  all  the  methods  discussed  here  and 
intercompare  the  results.  If  the  estimates  agree, 
they  can  be  recorded  with  great  confidence.  On  the 
other  hand,  the  difference  between  estimated  Rb  val¬ 
ues  from  R(z)  and  RL{z)  might  serve  as  a  crude  mea¬ 
sure  of  the  degree  to  which  the  bottom  behaves  as  a 
Lambertian  surface.  However,  for  purposes  of  mod¬ 
eling  the  light  field,  Mobley22  has  demonstrated  that 


measuring  the  magnitude  of  the  effective  bottom  al¬ 
bedo  is  far  more  important  than  obtaining  the  de¬ 
tailed  angular  pattern  of  the  bottom  bidirectional 
reflectance  function. 

Appendix  A:  Derivation  of  Eq.  (10) 

If  the  optical  properties  of  water  are  spatially  uni¬ 
form  and  there  are  no  internal  sources,  when  z  is  at 
least  one  optical  depth  away  from  a  boundary  the 
upward  radiance  and  downward  irradiance  can  be 
expressed  as  summations  of  eigenmodes13: 

Lu(z)  =  L(z,  -1) 

J 

=  2  [Cd^cKu,  -l)exp {-cz/vj) 
j=i 

+  C(-Vj)$(-Vj,  -l)exp(c2/u;)],  (Al) 

J 

Ed{z )  =  2  [C(U ) £ ! (o/ ) exp ( - C2 / U ) 

j=i 

+  Ci-Vjjg^-Vjiexpicz/Vj)],  (A2) 

where  C(±Vj)  are  expansion  coefficients,  v;  are  the  J 
eigenvalues  of  Eq.  (1)  corresponding  to  the  eigenfunc¬ 
tions23  4>(±Uy,  p),  and 


gj(±Vi)  =  f  4>(±u,  p)R/(p)d|x  (A3) 

Jo 

for  Legendre  polynomial  P7(p).  Far  from  the  bound¬ 
aries,  Ed(z)  and  Lu(z)  can  be  approximated  by  retain¬ 
ing  only  the  asymptotic  decreasing  eigenmode 
corresponding  to  the  largest  eigenvalue  vA  [i.e.,  C(vj) 
=  0  fory  >  1  and  C(—Vj)  =  0  for  all  y],  and  since  RL(z) 
=  tt Lu(z)/Ed[z)  the  asymptotic  value  of  RL(z)  is  R,L 
=  TTc^idi,  -1  )/g1{u1).  At  depths  far  from  the  surface 
but  where  some  influence  of  the  bottom  is  present, 
Ed(z)  and  LJz)  are  better  approximated  by  also  in¬ 
cluding  the  eigenmode  corresponding  to  the  largest 
negative  eigenvalue,  ~v1,  so  that 


Lu(z)  »  C(ui)<Koi,  —  l)exp(— cz/ui) 

+  C(-u1)4»(-i;1,  -l)exp(c2/u1),  (A4) 

Ed(z ) »  CivJg^vJexpi-cz/vJ 

+  Ci-vjg^-vjexpicz/vi).  (A5) 

After  forming  the  ratio  RL(z)  =  tt Lu(z)/Ed(z),  divid¬ 
ing  through  by  C(u1)^1(o1)exp(-cz/r;1),  letting  r  = 
C(— UiVCK),  and  recognizing  that13  Rj L  =  4»( — 
l)/£ifai),  we  find  that 


_  Rj1  +  urM-ib,  ~  l)exp(2cz/u1)/g1(u1) 
1  +  rR.,  exp  (2 Icz/iq) 

Subtraction  of  R.J  from  Eq.  (A6)  gives 


(A6) 


Rl(z)  -  Rj-  = 

r  exp(2czAq)[TT(K— u1;  ~l)/gi(vi)  -  RjR J 

1  +  rR,  exp(2cz/c,) 


(A7) 
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A  rearrangement  of  Eq.  (A7)  is 

[. RL(z )  -  E^exp^ 2cz/u1)[l  +  rRx  exp(2cz/o1)] 

=  r[m(>(— i>n  -l)/gi(i>i)  ~  RJ'R*].  (A8) 

Since  the  right-hand  side  of  Eq.  (A8)  is  independent  of 
z,  the  left-hand  side  at  arbitrary  depth  2  equals  that 
at  the  bottom.  Therefore 


RL(z)  =  R.,L  +  [RL(zb)  -  i?JJ]exp[-2c(zt  -  z)/iq] 
1  +  rR „  exp(2cz6/u1) 


X 


1  +  rR  ,  exp(2cz/u1) 


(A9) 


If  [ rR *,  exp(2czfe/r;1)]  <K  1  or  if  (zb  -  z)  is  small,  Eq. 
(A9)  reduces  to 


RL{z)  =  Rj*  +  [RL(zb)  -  i?0=t]exp[-2c(26  -  z)/v  1], 

(A10) 

which  is  analogous  to  our  equation  for  R(z)  derived  in 
a  similar  manner.17  Rearrangement  of  Eq.  (A10) 
gives  Eq.  (10). 
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